One-dimensional Bose chemistry: effects of non-integrability 
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Three-body collisions of ultracold identical Bose atoms under tight cylindrical confinement are 
analyzed. A Feshbach resonance in two-body collisions is described by a two-channel zero-range 
interaction. Elimination of the closed channel in the three-body problem reduces the interaction 
to a one-channel zero-range one with an energy dependent strength. The related problem with an 
energy independent strength (the Lieb-Liniger-McGuire model) has an exact solution and forbids 
all chemical processes, such as three-atom association and diatom dissociation, as well as reflection 
in atom-diatom collisions. The resonant case is analyzed by a numerical solution of the Faddeev- 
Lovelace equations. The results demonstrate that as the internal symmetry of the Lieb-Liniger- 
McGuire model is lifted, the reflection and chemical reactions become allowed and may be observed 
in experiments. 
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Ultracold atoms under tight cylindrical confinement 
can reach the "single-mode" , or quasi-one-dimensional 
(ID) regime, where only the ground state of trans- 
verse motion is significantly populated. Quasi-ID atomic 
quantum gases have been realized recently in elongated 
atomic traps (see HI 0, 0, Q), two-dimensional optical 
lattices (see [a, IE 13), atomic waveguides (see @ and 
references therein), and atomic integrated optics devices 
(see 0, ^3 and references therein). These systems at- 
tract increasing attention due to their possible applica- 
tions to atomic interferometry, quantum measurement, 
and quantum computation. 

Properties of quasi-lD systems are frequently analyzed 
using the Lieb-Liniger-McGuire (LLMG) model [lg, LL3] 
with S function atom-atom interaction U a S (z), where z 
is the interatomic distance. The interaction strength U a 
has been related in Ref. 01 (being denoted as gm) to 
the three-dimensional (3D) elastic scattering length a^o 
and the waveguide transverse frequency u>± . The LLMG 
model has an exact Bethe-ansatz solution 0, 0] ex- 
pressed as the superposition of plane waves with all pos- 
sible permutations of the asymptotic momenta pj, one 
per each atom. Therefore, the atoms can exchange their 
momenta, but the asymptotic momentum set {pj} re- 
mains unchanged. 

If the interatomic interaction is attractive (U a < 0), 
the atoms can form a two-body (2B) bound state (di- 
atom), described by complex asymptotic momenta p±in, 
where k = — y£/ a serves as a reciprocal diatom size and 
m is the atomic mass. However, a three-body (3B) as- 
sociation is forbidden in the LLMG model, as it cor- 
responds to a change of two real momenta in the set 
{pi,P2,P3} for a three-atom system by complex ones in 
the set {— p,p/2 + in, p/2 — in}. The dissociation and 
reflection in atom-diatom collisions are forbidden by the 



same reason (the last process corresponds to a transition 
{— p,p/2 + in, p/2 — in} — > {p, —p/2 + in, —p/2 — in}). 
Therefore, all chemical processes, i. e., atom-diatom 
transitions are forbidden in the LLMG model. 

The integrability of the LLMG model results from its 
high internal symmetry. This symmetry can be violated 
in real physical systems, e. g., when Feshbach resonance 
is used to tune the elastic scattering length. The effect 
of Feshbach resonance appears in 2B scattering when the 
collision energy of a pair of atoms in an open channel is 
close to the energy of a bound state (resonant molecule) 
in a closed channel (see Ref. ^4|). Two-body Feshbach 
resonance collisions in atomic waveguides have been ana- 
lyzed in Ref. [15J . The results demonstrate that when the 
collision energy is low enough compared to the trap fre- 
quency, the problem can be treated as essentially ID, but 
with an energy-dependent zero-range interaction. The 
present work demonstrates that a Feshbach resonance 
leads to non-integrability of few-body problems, allow- 
ing for chemical processes (association and dissociation) , 
as well as atom-diatom reflection. 

The approach used in Ref. ^{| for two-atom scatter- 
ing leads in a 3B problem to a system of two coupled 
equations for the wavefunctions in the momentum repre- 
sentation, ipo (gi, 52 7 93) and ip\ (91, 9 m ), associated with 
the three-atom and atom-molecule channels, respectively. 
Using a system of units in which % = 1, the equations are 



Epo (^1,92,93) = -7T- y^Vo (91,92,93) 
2m ' J 

3 

+ ^ U « E / d3 l' S - 9,) 5 (Q - Q 1 ) tpo (9i, 9 2 , 9 3 ) 



2 



Eipi (<?i,«?m) 



(37T) 

2m 

1/2 



-1/2 



9m 

4m 



^2<Pi (ijiQ-Qj) (!) 

.7 = 1 



ID 



<Pi (?l,?n 



for the symmetric transition amplitude X(p,po). Here 
the 2B scattering amplitude Ti£> (9), corresponding to 
the momentum-dependent U c g (q), has the form (see Ref. 

0) 
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Here qj are the momenta of the atoms, q m is the mo- 
mentum of the resonant molecule, and Q = 91 + 92 + 9.3 
is the center-of-mass momentum. The parameters of the 
ID problem — the atom- atom interaction strength U a , 
the atom-molecule coupling strength <?, and the detun- 
ing between the atomic and molecular states Did — can 
be expressed in terms of 3D resonance scattering pa- 
rameters (see Ref. (lif). The indistinguishability of the 
bosonic atoms leads to the symmetry of the wavefunction 
ipo (§1,92)93) over permutation of the atomic momenta. 
The energy E is counted from the threshold of the three- 
atom channel. 

Equation neglects the interaction associated with 
collisions between the resonant molecule and the third 
atom. It allows a simple elimination of the atom-molecule 
channel function ipi. As a result, the three-atom channel 
wavefunction obeys the equation 
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Eq uations similar to Eq. (J5J have been used in Refs. 
|20j for non-resonant ID problems. 

The poles of T 1D (q) at q = ±2y/m(E - E b )/3 cor- 
respond to 2B bound states (diatoms) with the binding 
energy E b — n 2 /m. The reciprocal diatom size k here 
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is a solution of the cubic equation (see 
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involving an effective atom-atom interaction strength 



U eS (9) = U a + 



E + iO-D 1D -3q 2 / (4m)' 



(4) 



This strength appears to be the same function of the col- 
lision momentum as in the 2B problem (see Ref. [lif). 
but here the collision energy of two atoms is expressed 
in terms of the total energy E and the third atom mo- 
mentum 9 due to momentum and energy conservation 
(letting Q = in the center-of-mass system used here 
and below). 

For a momentum dependent U e g the internal symme- 
try of the LLMG model is broken, and a resonant case 
does not allow a Bethe-ansatz solution. This case is 
analyzed here by a numerical solution of the Faddeev- 



These states are superpositions of the closed and open 
channels. The probability to find the diatom in the open 
channel is Wo = (k 2 + mDm) / (3k 2 + mU a n + mD\D) ■ 
This probability varies from 1 far off resonance to 3 in 
resonance for U a = 0. 

In the numerical solution of Eq. |5| the transition am- 
plitude is expressed in terms of odd and even amplitudes 
X± (p,Po) = X (p,p ) ± X (-p,po), p > 0, that satisfy 
uncoupled equations. The contributions of the poles of 
Tid (9) are separated from the Cauchy principal-value 
integrals as in Ref. [2(|. At E > 0, when chemical pro- 
cesses are allowed, the function Z (q,p) has singularities 
on the real axis too. Unlike the 3D case with logarithmic 
singularities, in the ID case the singularities are simple 
poles and are treated in a similar way as the poles of 

TlD(q)- 

The reflection, transmission, and dissociation proba- 
bilities in the collision of an atom with momentum po 
and a diatom with momentum — pq can be expressed in 
terms of the amplitude X (p,po) as, respectively, 
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Lovelace equations. Usin 
duction method (see Ref. 
equation 



the conventional Faddeev re- 



18]), one obtains a ID integral 
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Here the factor Wq describes the open-channel fraction 
in the incoming and outgoing waves. The results can 
be conventionally expressed in terms of dimensionless 
parameters: the non-resonant interaction strength u = 
m i/3|^|-2/3jj Qj the collision energy e 2 = 3pg/ (4m£> ), 
and the detuning b = Dm/ Do, where an energy scale 
Do = m 1 / 3 ^! 4 / 3 is used. 

Consider an association of three atoms with momenta 
Pi , P2, and P3 (pi + P2 + P3 = in the center-of-mass 
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FIG. 1: Probabilities of reflection (solid line) and dissociation 
(dashed line) in atom-diatom collisions calculated for b = 
and u — as functions of the dimensionless collision energy 
£2- The scaled rate coefficient of 3B association K^m/h, cal- 
culated as a function of the three-atom energy £3, is presented 
by the dot-dashed (it = 0, b = 0, 9 — ir/6) and dot-dot-dot- 
dashed (u = 1.6, b = 1, 6 = 7r/10) lines. 



±p* are determined by energy conservation, as 

E=±(pl+ P l+pl)=^-pt-E b . (9) 
The association rate coefficient is given by 
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Except for the dimensionless three-atom collision energy 
£3 = E/Dq, the rate coefficient depends on the energy 
distribution between the three atoms. This distribution 
can be parameterized by an angle < 9 < 7r/6, such that 
pi = y/3mE/4cos8, p 2 , 3 = ^SmE/Acos (9 ± 2tt/3). At 
9 = one of the momenta attains the maximal possible 
value pi = y/SmE/A, while 7^,3 = — Pi/2. A symmetric 
configuration (pi = —pi = \j mE, p% = 0) corresponds to 
9 = 7r/6. The association rate coefficient is plotted below 
in units of H/m. 

The calculated probabilities and rate coefficient are 
presented in Fig. as functions of the corresponding di- 
mensionless collision energies. The results demonstrate 
total reflection for slow collisions. Similar reflection 
behavior has been obtained in Ref. [13] for a case of 
non-identical atoms with energy-independent interaction 
strengths. The dissociation appears above the threshold, 
where 3pg/ (4m) > Ef, (or £2 > 1 for the conditions in 
Fig. ^) . The 3B association rate attains its maximum at 
a collision energy £3 ranging between 2 and 4 (depending 
on u, b, and 9), and decreases to zero for slow collisions. 
These processes vanish in fast collisions, when the inter- 
actions become small perturbations. The dependence of 



FIG. 2: The scaled rate coefficient of 3B association K^m/h 
calculated as a function of 6 for £3 = 2.4, u = 0, b = (solid 
line), £3 = 4.5, 6 = 1, u = 1.6 (dashed line), and £3 = 4.5, 
b = 1, u = —2 (dot-dashed line). 




FIG. 3: Probabilities of reflection (solid line) and dissociation 
(dashed line) in atom-diatom collisions calculated for the fixed 
dimensionless collision energy £2 = 3.6 and the detuning b — 
as functions of the dimensionless detuning b. The scaled rate 
coefficient of 3B association K^m/h is presented by the dot- 
dashed (u = 0, £3 = 2.4, 8 — 7r/6) and dot-dot-dot-dashed 
(u = 1.6, £3 = 4.5, 6 = tt/10) lines. 



the association rate on the energy distribution between 
the atoms behaves differently for various u and b values 
(see Fig. 0). No association occurs at 9 = when two 
atoms have equal momenta. 

Given a fixed collision energy, the reflection, dissoci- 
ation, and association persist over a detuning interval 
(see Fig. |3J). The oscillating behavior of the probabili- 
ties results from interference of the odd and even ampli- 
tudes X± (p,po), which are analogs of partial waves in 
3D scattering. However, unlike the 3D case, where only 
the s-wave contributes to low-energy scattering due to 
a centrifugal barrier, in the ID case both odd and even 
amplitudes have to be taken into account at all energies. 
All chemical processes, as well as reflection, vanish at 
large detunings, when the energy dependence of the in- 
teraction strength U g becomes negligible [see Eq. (0}], 
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FIG. 4: Probabilities of reflection (solid line) and dissoci- 
ation (dashed line) in atom-diatom collisions calculated for 
fixed dimensionless collision energy €2 = 3.6 and detuning 
b — as functions of the dimensionless non-resonant interac- 
tion strength u. The scaled rate coefficient of 3B association 
Ksm/h is presented by the dot-dashed (b = 0, 63 = 2.4, 
9 = tt/6) and dot-dot-dot-dashed (6 = 1, e 3 = 4.5, 6 = it/10) 
lines. 



restoring integrability of LLMG model. For a repulsive 
non-resonant interaction (U a > 0) the bound state exists 
while Did < 2\g\ 2 /U a [see Eq. 0], or while b < 2/u. As 
b — > 2/u the bound state becomes more shallow, leading 
to a drastic increase of the association rate. The depen- 
dence on the non-resonant interaction strength (see Fig. 
|3J) demonstrates a similar behavior : a drastic increase of 
the association rate near the threshold at u = 2/6, and a 
decrease of all the rates at large absolute values of u due 
to the restoration of integrability. 

Thus, a substantial association rate can be observed 
for collision energies ranging over several units Do, which 
must lie below the transverse frequency in order to 
conserve the quasi-lD behavior of the system. This con- 
dition is satisfied for relatively weak resonances /xA < 
w±a±/a3D and small detunings /i (B — Bq — A) — lo± < 
y // IJ]jlAa±Ja^ , when the relations between the ID and 
3D scattering parameters of Ref. ^5| are applicable. 
These relations lead to expressions for the energy unit, 



1/2 

its resonant value Bo, a± = [2/ (muj±)] is the trans- 
verse harmonic oscillator length, and C w 1.4603 (sec 
Ref. ^2). For example, in a waveguide with the trans- 
verse frequency u>± — 50 x 27r KHz, collisions remain 
quasi-lD for the collision energy less then 2.4//K. In this 
case D Q = 0.27/xK, u = 0.22, dB/db = 1.1 mG, and 
%/m — 2.6 x 10 -5 cm 2 /s for the Na resonance at 853 
G; D = 0.48^K, u = 0.53, dB/db = 0.52 mG, and 
h/m — 6.9 x 10~ 6 cm 2 /s for the 87 Rb resonance at 685 
G; and Do = 0.35^K, u = 1.7, dB/db = 9.3 mG, and 
Ti/m = 4.5 x 10~ 6 cm 2 /s for the Cs resonance at 20 G. 

The high internal symmetry of the LLMG model can 
be also violated due to virtual transverse mode excitation 
even if the interaction of free atoms is non-resonant. This 
effect has been considered in Ref. 22] for 3B collisions in 
atomic waveguides and in Ref. |23j for atom-soliton col- 
lisions. It can be also analyzed using the present ap- 
proach, since a confinement-induced resonance can be 
represented as a kind of Feshbach resonance (see Ref. 
|12j). The scattering amplitude in an atomic waveguide 
see Eq. (6.9) in Ref. can be approximated by Eq. 
(|rj|) with the proper choice of the parameters g, Did, and 
U a . However, now the problem contains the single energy 
scale, uj± , and the energy unit Dq will have the same or- 
der. As a result, the association and dissociation rates 
are negligibly small in the quasi-lD regime and the only 
evidence of nonintegrability remains the atom-diatom re- 
flection, in agreement with Refs. jU|2^|. A narrow Fesh- 
bach resonance can provide reacher physical phenomena. 
In addition to association, dissociation, and reflection the 
nonintegrability can lead to a change of the asymptotic 
momentum set in three-atom elastic collisions. 

In conclusion, the 3B problem with two-channel 2B 
interactions can be reduced to a single-channel case with 
an energy-dependent interaction strength. The energy 
dependence violates the high internal symmetry of the 
LLMG model, allowing chemical reactions. Reflection 
becomes the dominant output channel of atom-diatom 
low-energy collisions, while dissociation appears above a 
threshold. The 3B association should be observable in 
experiments. 
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and the dimensionless parameters 
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Here A is the phenomenological resonance strength, \x 
is the difference between the magnetic momenta of an 
atomic pair in the open and closed channels, B — Bq 
is the detuning of the external magnetic field B from 
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